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Decay of a metastable state activated by non-Gaussian noise: 
A critical review of the generalized Kramers problem 
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We review the problem of the noise activated escape from a metastable state in the presence of non- 
Gaussian noise, and present connections between various theoretical approaches. We also respond 
to criticism of our work by Tomas Novotny [arXiv;0807.0387] concerning the weak damping limit. 
The discrepancy between our results is linked to uncontrolled approximations made by Novotny. 
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Introduction. — The problem of noise activated escape 
from a metastable state has been formulated and solved 
by Kramers more than 50 years ago [1]. This important 
problem has many appUcations across different fields. It 
has been addressed many times in the literature, and re- 
viewed in many textbooks, e.g. in [2]. By now, this prob- 
lem is well understood on the level of Gaussian noise. Re- 
cently, new interest in Kramers' problem has been stim- 
ulated by experimental progress in the field of electronic 
on-chip detectors, which allows the measurement of small 
non-Gaussian effects in electronic noise. Kramers' prob- 
lem was reconsidered beyond Gaussian noise and formally 
solved in Refs. [3-5] for the case of strong damping and in 
Ref. [6] for the case of weak damping, where the number 
of variables is reduced. This problem has been applied 
in Refs. [6-8] to the stochastic dynamics of the Joseph- 
son junction (JJ) threshold detector, where the escape 
rate from the supercurrent state of the JJ provides the 
information about the current noise. 

It has been shown that the third cumulant of cur- 
rent noise, which describes a weak non-Gaussianity of 
a Markovian process, contributes to the asymmetry of 
the escape rate with respect to the current bias. Very 
recently, this phenomenon has become a point of con- 
troversy. In the paper "Josephson junctions as thresh- 
old detectors of the full counting statistics: open is- 
sues" [9] Tomas Novotny claims that his new theory 
of noise-activated escape makes a prediction concerning 
non-Gaussian effects which contradicts previous theoret- 
ical results [6, 8]. The disagreement only concerns a 10% 
difference in a detector parameter's coefficient, and there- 
fore would not otherwise be of great concern. However, 
Novotny goes on to comment that "...it is relevant from 
a purely conceptual point of view which one is actually 
correct since it should help with the identification of pos- 
sible misconceptions hidden in the failed approach (es)." 
The purpose of the present short paper is to clarify and 
resolve this controversy and to explain the connection 
between different theoretical schemes. 

Here we show that the Novotny theory represents a 
particular case of the stochastic path integral (SPI) for- 
mulation introduced in the author's earlier paper [6] and 



relies on an alternative method of solving the resulting 
Hamilton- Jacobi equation for the classical action. How- 
ever, the presented solution is not correct, which explains 
the discrepancy with previous theoretical results. The 
mistake is not simple, but concerns the subtle nature of 
the weak-damping limit and an uncontrolled approxima- 
tion made by Novotny. Rather than simply point to the 
mistake, we take this opportunity to review the physics, 
and show how the approaches used in several papers are 
related. 

We start with the brief formulation of the problem. 
Consider a particle which moves in a metastable poten- 
tial V{q) around the local minimum at point q = qo (see 
Fig. 1). Such a particle may represent a system in a 
state described by the collective variable q. The conser- 
vative classical motion is generated by the Hamiltonian 
H{p, q) = p^ /2m + V{q). In addition, if the particle in- 
teracts with the environment, this leads to damping, so 
that the particle relaxes to the point at the minimum of 
the potential, q — qo. At the same time, the environmen- 
tal noise activates the motion, so the particle may escape 
from the metastable state via the point q = qi. We fur- 
ther assume that this complex behavior may be described 
by the set of Hamilton-Langevin (HL) equations, 

q = dH{p,q)/dp + I„ p=~dH{p,q)/dq + Ip, (1) 

where the currents Iq and Ip are the sources of noise, 
which also describe damping. The problem is to find the 
rate of escape. 

Formal solution of Kramers ' problem. — If the currents 
are fast variables, so that they fluctuate on the time 
scale To, which is much shorter than the the characteristic 
time scale T of the deterministic motion of the particle, 
then the HL equations can be solved by the method of 
the stochastic path integral, introduced for non-Gaussian 
noise in Refs. [10, 11]. According to this method, the evo- 
lution of the particle is described by the path integral 

P= fvA fvRexp{S), (2) 

where the action S is given in the explicitly canonically 




FIG. 1: A metastable system with the collective coordinate 
q can be viewed as a "particle" which moves in a potential, 
shown by the full line, with a local minimum at the point 
q = qo and a maximum at the point q — qi. The dotted line 
schematically shows the most likely trajectory which leads 
to the escape of the system from the metastable state. In 
the weak damping limit the motion is quasi-periodic, and the 
number of winding increases with the quality factor Q of the 
system. 



invariant form as 

S= f dt'[-A-R. + A- {R, H} + n{k, R)]. (3) 

Here R = (p, q) and A = (Ap, A^) are the sets of physical 
and canonically conjugated auxiliary variables, respec- 
tively, {. . .} denotes the Poisson bracket with respect 
to p and q, and the function ?i(Ap,Ag) is the gener- 
ator of the cumulants of the Markovian noise sources: 
((/^/™)) = dl^d-^niO, R) and n{0, R) = 0. By fixing R 
in the final state of (3) we obtain the probability distri- 
bution P(R), while fixing the final A variables turns the 
SPI into the moment generating function P(A). That is 
why the auxiliary variables A are also called "counting" 
variables in the field of full counting statistics [12]. 

The separation of time scales T ':S> T{) leads to a small 
parameter tq/T which plays the role analogous to h in the 
Feynman path integral and thus allows the saddle-point 
evaluation of the SPI. This naturally leads to Hamilton's 
equations of motion in the extended space [11] 

R = dK/dk, A = -dK/dH, (4) 

generated by the new Hamiltonian 

isr(A,R)EEA-{R,i/} + W(A,R), (5) 

which has to be solved for some initial condition R(to) = 
Ro and A(to) = Ao and substituted back to the action 
to give the evolution operator Pit) = exp[5(t)]. 

In the context of Kramers' problem there always exists 
a trivial solution A = and R = {R, iJ} + (I), where 
I = (/p, Iq), for the "average" dynamics in physical space 
with a null Hamiltonian and action S* = -fC = 0, giving 
the proper normalization of the distribution P. This so- 
lution describes the relaxation to the local minimum of 
the potential at g = go (resting state), and away from 



the top of the barrier for q > qi (the running state). We 
are looking for a nontrivial solution with non-zero action 
which connects these two states and thus leads to the es- 
cape from the local minimum (it is shown in Fig. 1 by a 
dotted line). Therefore, the initial condition for this so- 
lution is A(to) ~ 0, p{to) = 0, and (/(to) = <Zo- Moreover, 
since the full Hamiltonian is an integral of motion, and 
for the initial state oi K ^ 0, we have 



K[A{t),Ii{t)] =0 



(6) 



along the trajectory of interest. Therefore, the classical 
action simplifies and we find 
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dtA{t) ■ R(t) ^ - / A(R) • dR. (7) 



Rn 



Then, up to a prefactor [13], the rate of the noise ac- 
tivated escape from the localized state is given by T oc 
exp(S'). 

Although the original problem is stochastic, the ad- 
vantage of the SPI method [11] is that it reduces the 
problem to solving deterministic equations of motion in 
the extended phase space. Therefore, Eqs. (4-7) provide 
a unique solution of Kramers' problem for an arbitrary 
underlying conservative dynanrics generated by H{p, q) 
and for arbitrary Markovian noise generated by ?i(A, R). 
However, these equations have a different equivalent rep- 
resentation, which may have triggered the confusion in 
the paper [9]. In the next part of the paper we clarify 
the connection between these two representations. 

Hamilton- J acobi equation. — Following the textbook of 
Landau and Lifshitz [14], we note that according to Eq. 
(7) the classical action as a function of coordinates sat- 
isfies A = —dS/dR. We will use the symbol S to repre- 
sent both Hamilton's characteristic function (referred to 
as action), and (when evaluated between the time limits 
on the escape trajectory) the exponential contribution to 
the activation rate, as in (7). Following this observation, 
equation (6) can be written in the form of a differential 
equation for the action. 



K{-dS/dIl,Il) = 0, 



(8) 



which is nothing but the Hamilton-Jacobi (HJ) equation 
in the stationary case. A solution of this equation in 
some region R e f2 requires fixing a boundary condition 
at the edge dfl, which can be conveniently written as 
dS/dRlon — — Ao. There is nothing special about the 
direct solution of the equation (8), if the correct physical 
boundary conditions are satisfied. However, one should 
remember that such a solution is unique and should co- 
incide with the one that follows from the equations of 
motion. Indeed, one can equally consider the bound- 
ary condition for dS/dR as the initial condition for the 
canonical momentum A (to) = Ao and for the coordinate 
R(to) G dV,. Then the function K generates the Hamil- 
tonian dynamics with the canonical momentum A(R) 



which coincides with the vector field —dS/dR in the re- 
gion O. This procedure can be viewed as a method of 
characteristics for solving the equation (8). 

The way Novotny [9] arrives at the HJ equation de- 
serves a separate consideration. Here we present the most 
general variant of it. Let us consider a specific noise pro- 
cess, one which is generated by rare random transitions 
from the state R to the state R' with the rate r(R', R). 
Since the transitions are rare (more rigorously, Ttq <C 1), 
the probability distribution function P(R, t) satisfies the 
following master equation [2] 

P = dP/dR ■ R + dP/dt = {R, H} ■ dP/dH 

+ /"dR'[r(R,R')P(R')-r(R',R)P(R)], (9) 

where the first term takes into account the determinis- 
tic part of the dynamics. We follow then the standard 
Kramers-Moyal forward expansion [2] to derive the equiv- 
alent Fokker-Planck equation. The result is 

p=:{R,i/}-ap/aR + H(-aR,R)P(R), (lo) 

where the function Ti is defined as follows 

W(A,R) = /dRir(Ri,R) U^i-^)^ - ll . (11) 



In general, the order in which the derivative (9r, is taken 
in (10) is important. However, accounting further for 
the separation of time scales parameter T/tq ^ 1, and 
making the semiclassical approximation P ~ cxp(S'), we 
apply the derivatives in (10) directly to the action. In 
this way in the stationary case P = we again obtain 
the HJ equation (8) with K given by (5) and the function 
7i given by Eq. (11). The paper [9] then recognizes this 
function as the generator of the Poissonian process and 
generalizes the result for an arbitrary noise. This step, 
although it gives the correct HJ equation, is not rigorous. 
In particular, this step cannot be directly applied to the 
Fokker-Planck equation (10), which simply does not exist 
for the case of general noise (or, if formally written, does 
not correspond to any physical reality) . This is because a 
particular form of the FP equation fixes the quantization 
procedure via the operator ordering in (10), while the 
function Ti, does not provide enough information for this. 

Alternatively, on the SPI level the FP equation may 
be obtained from (2) in a standard way by expanding 
the "one-step" propagator P(R, R'; Ai) with respect to 
Ai. This step formally requires AtTi to be small. On 
the other hand, the SPI is constructed with the separa- 
tion of time scales requirement, T ^ Af 3> tq, which 
inevitably leads to the restriction Tirg <C 1 . In turn, this 
implies that the noise is generated by a weak Poissonian 
process, Pro <C 1, with the generator (11). However the 
SPI method itself does not rely on writing the FP equa- 
tion as an intermediate step, and therefore it is free from 
this difficulty and provides a rigorous and general way to 



describe the stochastic system on the semiclassical level. 
This fact is thoroughly explained in Ref. [11]. 

Energy diffusion. — To conclude the formal discussion 
of the problem, we consider the weak damping limit. In 
this limit the system experiences a quasi-periodic motion 
(see Fig. 1). The energy of the system weakly changes 
over the period making small steps back and forth. This 
process then can be viewed as diffusion of the energy, 
so that the action depends only on the energy variable. 
For the rigorous derivation of the action via the canoni- 
cal transformation of coordinates we refer the reader to 
our paper [6], and present here a simplified physical ar- 
gument. 

The energy balance equation 



E = qlp- piq 



(12) 



follows directly from the HL equations (1) and takes the 
form of a Langevin equation. If Tq 3> tq, then variables 
q and p change slowly and can be considered as "effec- 
tive charges" . Therefore, the noise source on the right 
hand side of the equation (12) is described by the cumu- 
lant generating function H^XeQ, —^ep)- Diffusion of the 
energy is then described by the SPI with the action 



dt'[-XEE + H{XEq,-XEP)]- 



(13) 



To leading order in weak damping we can replace the 
generator Ti. in Eq. (13) with its average over the period 
T of oscillations 

{n)E = T-i i dtH{XEq, -Xep), (14) 

evaluated for fixed A^; and E. Corrections in damping 
will be found by taking into account slow energy dissipa- 
tion E ~ dx^H and A^; = — SbH, while averaging over 
the period T. This gives the action 



S 



dt'{-XEE+{H)E)- 



(15) 



The simplification here arises from the fact that in con- 
trast to the original action (3), the Hamiltonian K — 
{'H)e in the action (15) depends on one variable only: 
the energy E. Therefore the condition K ~ becomes 
an algebraic equation for the instanton line Xe — Xm (E) , 
which connects the bottom and the top of the potential 
barrier, see Fig. 1. Thus we arrive at the following result 
[6] 



S = - XindE, {H{Xinq, -Ai„p))_E = 0, 



(16) 



which formally solves the Kramers' problem for an arbi- 
trary Markovian noise in the weak damping limit. This 
result will be compared with the alternative way of ob- 
taining the weak damping limit. 



Josephson junction threshold detector. — The idea to 
use a Josephson junction in the metastable supecurrent 
state as a detector of noise belongs to Pekola [15] and 
to Tobiska and Nazarov [3]. Later on, the problem was 
reconsidered on a different level in the papers [6-9] . Leav- 
ing the experimental realization [16-18] of the idea aside, 
we focus solely on the theoretical part of the problem. 
For the physics of the JJ detector we refer the reader to 
our paper [6]. Here we just mention that after rescaling 
the physical variables [9] the dynamics of the detector is 
described by the HL equations (1) for a "particle" with 
mass m ~ I and with coordinate q (being a supercon- 
ducting phase), which moves in the periodic potential 
biased by the supercurrent J: 

H{p,q)=py2 + V{q), V{q)^-COs{q)-Jq. (17) 

The dissipative part of the current localizes the JJ in 
the supercurrent state in one of the local minimums of the 
potential, e.g. at q — qa, while the current noise activates 
the JJ and leads to the escape to the running dissipative 
state. This effect enters via the momentum source Ip, 
while the coordinate source Iq ~ 0. Truncating the noise 
generating function at the third cumulant, we write 



'^ =" -F;i-P>^P + ^p) + C3>^1 



(18) 



where Q is the quality factor of the JJ oscillator in the 
supercurrent state, and the constant C3 is proportional to 
the third current cumulant, ((/p)) = d'l Ti — 6C3, which 
characterizes a weak non-Gaussianity of noise. In all the 
papers on the subject C3 is considered to be a small ex- 
pansion parameter. 

Finally, the HJ equation (8) acquires the following 
form: 



dS 



V'iq) 



dS_ 
dp 



P_dS_ 
Qdp 



1 



dS_ 
dp 



= C3 



dS_ 
dp 



(19) 
where we deliberately neglect important physical effects 
in order to arrive at the HJ equation in the form pre- 
sented in Ref. [9]. Indeed, in general the source Ip de- 
pends on the system's state, which leads to the "cascade 
corrections" [11] to the action. In the context of the J J 
physics, this is an effect of the measurement circuit back- 
action. Reference [6] fully accounts for these circuit ef- 
fects and comes to the conclusion that the corrections are 
small and can be neglected if the load circuit impedance 
is small. This, in fact, is the case in all existing experi- 
ments. 

Weak non-Gaussianity. — We are now in the position 
to critically compare the results of Refs. [6-9]. All of the 
papers implement a perturbation expansion in the small 
parameter C3 to leading order. We start with our origi- 
nal work [6], where the third cumulant contribution was 
found in the strong damping Q ^ and the weak damp- 
ing Q ^ 00 limits. We focus here on the weak damping 



limit and reduce the dynamics in the full space to the 
diffusion in the energy space, as described above. Using 
Eqs. (16) and (18), we solve the equation for the instan- 



ton line perturbatively in C3, Ajn = 
and arrive at the following result 



1~csQ{p^)e/{p^)i 



s - -Ay + C3S3 



Vi 



{P')i 



S3^Q dE ■ , 

Jvo {P)E 



(20) 



where AV = Vi — Vq is the potential barrier, i.e. the 
first term is just Kramers' result for the Gaussian noise. 
Note that the second term contains a large prefactor Q — > 
00 which, on the other hand, is multiplied by the small 
integral, because {p^)e = to leading order in damping. 
This can be explained by the fact that for large Q the 
escape trajectory (shown in Fig. 1) makes small energy 
steps and large number of windings. In the appendix of 
Ref. [6] we take weak damping into account, E = p^ /Q 
(as explained above), to find the next order correction to 
{p^)e, and arrive at 



Vi 

83 = 2 I dE 

ha 



{W-<i)p')l 

{p')e 



(21) 



where q' is one of the two turning points of the conserva- 
tive motion, which is closest to the point qi . Thus we find 
that the quality factor Q cancels from the final result. 

Grabert's paper [8], despite being long, can be summa- 
rized in just one paragraph. It utilizes the SPI method of 
Ref. [6] and directly solves Eqs. (4-7) for arbitrary damp- 
ing. On the Gaussian level the escape path that satisfies 
Hamilton's equations of motion and all the boundary con- 
ditions listed before Eq. (6) is described by the equations 

q = P, A, = y', p^-r+p/Q, Xp^p, (22) 

which is time reversed with respect to the corresponding 
relaxation path. Substituting these equations into the 
action (7), we obtain the Gaussian part of the action. 



SG(p,q) = -H{p,q) = -pV2 - V{q), 



(23) 



which when evaluated between t^ and ti along the es- 
cape trajectory gives Kramers' result Sq = — AV^. Next, 
considering the term Ti' = csA^ in (18) as a perturba- 
tion, we note that to leading order in C3, the variation 
of the trajectory does not contribute to the action (3), 
so 5 = ^G + / dtJi' . Therefore, using again Ap = p we 
obtain that the total action takes the form (20) with 



& = 



dt[p{t)]\ 



(24) 



to 



where the integral is taken over the escape path. This 
is the main result of the Ref. [8] , which generalizes the 
calculations of Ref. [6] to arbitrary damping Q. In the 
weak damping limit, one easily obtains Eq. (20). In- 
deed, we can write / dtp^ = / dt{p^)E = / dE{p^)E/E = 
Q j' dE (p^) E / (p^) E , where we have used E —p^/Q. 



Finally, the Refs. [7, 9] choose to solve the HJ equation 
(19) using different approximations. Here we reproduce a 
few steps of Ref . [9] . To leading order in the small param- 
eter C3 one can write S ~ Sq + £36*3. Substituting this 
expansion to the equation (19) and using the Kramers 
action (23), one arrives at the following equation: 



dq dp Q dp 



P 



(25) 



This equation can be solved exactly using the method 
of characteristics by choosing S3[p{t), q{t)] with p{t) and 
q{t) being solutions of the equations (22) and describing 
the escape path in Fig. 1. In this case, equation (25) 
may be rewritten as dS^/dt = p^ by using equations 
(22), which recovers the previous result (24). There- 
fore, the direct solution of the HJ equation gives ex- 
actly the same result as solving Hamilton's equations of 
motion, as expected. The mistake of Ref. [9] is that it 
treats the limit Q ^ 00 incorrectly by dropping the l/Q- 
term in the above equation, giving Novotony's equation 
pdqS^ — V'{q)dpSz = p^- This step would be equiva- 
lent to dropping the 1/Q term in the equations of mo- 
tion (22) in our approach. Indeed, the dropping of this 
term corresponds to an exactly constant energy, E — 0, 
rather than the physically correct slowly growing energy, 
E = p^ IQ. After this uncontrolled approximation, the 
Novotny equation may be integrated over the energy con- 
serving periodic trajectories, leading to ^3 = / dtq-f? . Af- 
ter integrating by parts and using p = —V one arrives at 
'S'3 = (91 — 9o)p^ + "^ J dqV (q)q, which is the expression 
(9) of the Ref. [9]. However, this solution corresponds 
to exactly energy-conserving trajectories, and therefore 
cannot be correct. This explains why Ref. [9] finds the 
parametrically correct contribution of the third cumu- 
lant, while the dimensionless function of the current bias 
Di{J) is wrong in the detailed dependence on J. 

To summarize, we have reviewed the generalized 



Kramers problem of the decay of a metastable state un- 
der the influence of the non-Gaussian noise. We com- 
pared different calculations and came to the conclusion 
that they utilize basically the same theoretical frame- 
work introduced earlier in the papers [10, 11]. However, 
Refs. [7, 9] make uncontrolled approximations in solving 
Hamilton- Jacobi equations and thus arrive at parametri- 
cally correct, but quantitatively wrong results. 
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